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ABSTRACT 

Insulin signaling in the CNS modulates satiety and glucose metabolism, but insulin target neurons are poorly defined. We have previously shown 
that ablation of insulin receptors (InsR) in Glut4-expressing tissues results in systemic abnormalities of insulin action. We propose that Glut4 
neurons constitute an insulin-sensitive neuronal subset. We determined their gene expression profiles using flow-sorted hypothalamic Glut4 
neurons. Gene ontology analyses demonstrated that Glut4 neurons are enriched in olfacto-sensory receptors, 1712 acetylcholine receptors, and 
pathways required for the acquisition of insulin sensitivity. Following genetic ablation of InsR, transcriptome profiling of Glut4 neurons 
demonstrated impairment of the insulin, peptide hormone, and cAIVIP signaling pathways, with a striking upregulation of anion homeostasis 
pathway. Accordingly, hypothalamic InsR-deficient Glut4 neurons showed reduced firing activity. The molecular signature of Glut4 neurons is 
consistent with a role for this neural population in the integration of olfacto-sensory cues with hormone signaling to regulate peripheral 
metabolism. 

© 2014 The Authors. Published by Elsevier GmbH. This is an open access article under the CC BY-NC-ND license (http;//creativecommons.org/licenses/by-nc-nd/3.0/). 
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1. INTRODUCTION 

Insulin acts through the insulin receptor (InsR) in the CNS to regulate 
food intake [1], counterregulatory response to hypoglycemia [2], and 
leptin sensitivity [3]. However, anatomic and functional mapping of 
insulin's CNS action is not nearly as detailed as those of other hor- 
mones [4]. This is partly due to the fact that InsR are thought to be 
widespread among neurons, not to mention glial, and non-neural cells 
[5]. Glucose utilization is thought to occur largely in an insulin- 
independent manner in the CNS [6]. Nonetheless, a sub-population 
of widely scattered neurons possesses the insulin-responsive 
glucose transporter, Glut4 [7]. While the function of this transporter 
in glucose metabolism continues to be a matter of speculation, it's a 
fact that removing InsR from Glut4 neurons has a profound detrimental 
effect on peripheral metabolism. Thus mice lacking InsR in Glut4- 
expressing tissues develop type 2 diabetes [7], highlighting the 
contribution of CNS insulin resistance to the pathophysiology of dia- 
betes [8]. 

In this study, we set out to investigate the nature of Glut4 neurons in 
the hypothalamus, a key regulatory site for appetite and glucose ho- 
meostasis [9]. As a first step toward characterizing the physiological 
function(s) of Glut4 neurons, we determined their gene expression 
profiles. To this end, we generated transgenic mice with chemically 
defined Glut4 neurons [7]. We flow-sorted Glut4 neurons and 
compared their gene expression profile with that of non-Glut4 neurons 
using a previously reported protocol [10]. Interestingly, we found 
pathways critical for sensory perception, receptor signaling, and 
response to nutrient are significantly enriched in the Glut4 neurons. 



consistent with the hypothesis that Glut4 neurons have critical roles in 
metabolic regulation. 

2. MATERIAL AND METHOD 

2.1. Mice 

Gt(Rosa)26Soi""^"^^'^-'^'^'""^'°>"'' mice were from the Jackson Labo- 
ratories. GIRKO mice were generated as previously described [7]. The 
Columbia University Animal Care and Utilization Committee approved 
all procedures. Normal chow diet (NCD) had 62.1% calories from 
carbohydrates, 24.6% from protein and 13.2% from fat (PicoLab ro- 
dent diet 20, 5053; Purina Mills). 

2.2. Flow cytometry, gene profiling, and quantification of neurons 
We dissociated mediobasal hypothalami from 3-week-old mice with 
transgenes {Glut4-Cre [line 535]; Rosa-Tomato\ InsR /ox allele) with 
papain dissociation kit (Worthington Biochemical). We gated live 
neurons to collect Rfp-positive neurons. We performed microarray for 
gene profiling using a procedure described previously [10]. We used 
Partek Genomic Suite and R for pathway analysis and heatmap gen- 
eration. We analyzed the FACS data using Flowjo software. We 
quantified immunofluorescence using Image J software. 

2.3. Electrophysiological studies 

For patch-clamp recording in acute slices, we sectioned 350-|im- 
thickn coronal brain sections in a cold sucrose-based cutting solution 
containing (in mM): 1 95 sucrose, 1 0 NaCI, 25 NaHCOs, 25 glucose, 2.5 
KCI, 1.25 NaH2P04, 2 Na pyruvate, 0.5 CaCb, and 7 MgCl2 (pH 7.3) 
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using a vibratome (Leica VTIOOOS). Slices were incubated in a sub- 
merged chamber at room temperature for at least 90 min before re- 
cordings, perfused with recording solution containing (in mM): 125 
NaCI, 2.5 KCI, 1.25 NaH2P04, 25 NaHCOs, 2 CaCl2, 1 IVIgCIa, 25 
glucose (pH 7.3) and continuously bubbled with 95% O2 and 5% CO2. 
We selected cells using an upright microscope fitted with fluorescence 
optics (Nikon) and patched under IR-DIC optics. We acquired patch- 
clamp recordings with a MultiClamp 700B amplifier and Clampex data 
acquisition software (Axon Instruments) at room temperature. We 
digitized data at 10 kHz with a 2 kHz low-pass filter. We fabricated a 
glass pipette of 5—8 IVIQ and filled it with electrode solution (in mM): 
130 K gluconate, 10 KCI, 5 HEPES, 5 IVIgCl2, 0.06 CaCl2, 0.6 EGTA, 2 
MgATP, 0.2 Na2GTP, and 20 Phosphocreatine. We adjusted osmolarity 
to 31 0 mOsm and pH to 7.3 with KOH. We assessed electrical access 
in the whole-cell recording mode using a small-voltage step pulse 
before and at the end of the current-clamp recordings. We discarded 
recordings in which the series resistance had been >30 MQ. Resting 
membrane potential and spontaneous action potentials were recorded 
in the current-clamp mode with no current injection. Detection and 
measurements of action potentials were performed using IVIiniAnalysis 
software (Synaptosoft). 

2.4. Statistical analyses 

We analyzed data with Student's f-test or 1 way or 2way ANOVA using 
GraphPad Prism and Partek Genomic Suite software. We used the 
customary threshold of p < 0.05 to declare statistical significance. 

3. RESULTS 

3.1. Glut4 neurons are enriched in olfacto-sensory receptors and 
hormone signaling pathways 

We took an unbiased approach to characterize Glut4 neurons. We 
generated Glut4-Cre: Rosa26-Tomato mice, in which Glut4 neurons 
were labeled by Cre-mediated activation of red fluorescent protein 
(Rfp). We used fluorescence-activated cell sorting (FACS) to collect 
Rfp-labeled Glut4 neurons and non-fluorescent control neurons from 
basal hypothalami of transgenic mice. We compared transcriptomes of 
Glut4 and control neurons by microarray analysis. Gene ontology (GO) 
enrichment analysis demonstrated that the Glut4 neuron transcriptome 
was enriched in olfacto-sensory systems and G-protein coupled re- 
ceptors. Other classes of enriched genes could be subsumed under 
these classifications. In addition, we found enrichment in genes 
required for nutrient sensing and glucose utilization, consistent with 
the notion that these neurons have metabolic control activity (Table 1). 
Indeed, GO ANOVA analysis indicated an overrepresentation of genes 
involved in the maturation of metabolically active cells (adipocytes). 
This class included transcriptional regulators such as Ppar-y, Tcf7l2, 
and others that could be subsumed under the classification of histone 
modification genes, which was also enriched in Glut4 neurons. GO 
pathway analysis revealed enrichment in genes required for 0-glycan 
biosynthesis and ribosome biogenesis (Table 1). 
In an attempt to assign Glut4 neurons to a specific class of neurons, we 
carried out a separate analysis of genes that are essential for neuronal 
activity, including neurotransmitter receptors and ion channels. A 
heatmap rendering shows that most neurotransmitter receptors have 
variable levels of enrichment in Glut4 neurons (Figure 1 A). We analyzed 
the data by grouping the major sub-classes of neurotransmitter re- 
ceptors (Figure IB— I). Overall, we found no enrichment of specific 
neurotransmitter subsets, indicating that Glut4 neurons likely include 
different neuronal subtypes (Figure IB— I). However, we found a sig- 
nificant enrichment of the M2 subtype muscarinic receptor, along with 
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Gene ontology enrichment 



Function 


Enriciiment 


p-Vaiue 


%Genes in group 




score 






represented 


Olfactoiy receptor acfi\flty 


17,3227 


3.00E-08 




4,07031 


Sensory perception of smell 


16.486 


6.93E-08 




3,96328 


Sensory perception of cherrrical stimulus 


16,3667 


2.14E-07 




3,76313 


Sensory perception 


11,6164 


9.03E-06 




3,20382 


G-protein couplerJ receptor activity 


10,9879 


1 .69E-05 




3,00678 


Neurological system process 


10,0768 


4.20E-05 




2,94638 


Transmembrane signaling receptor 


9,06199 


0,000115992 




2,74177 


activity 










G-protein coupier] receptor signaling 


8,72103 


0,000163119 




2,7376 


pathway 










Sign^ing receptor activity 


7,98683 


0,000339911 




2,62781 


Ceiluiar response to nutrient 


7,31929 


0,000662631 




27,2727 


Regulation of transcription by glucose 


5,66087 


0,00388407 




33,3333 


Gene ontology ANOVA 


GO description 


# of probes 


GO ID 


/^Vaiue (attribute) 


Fat ceii differentiation 




63 


45444 


0.00140471 


PfiospholipirJ homeostasis 




8 


55091 


0.00212119 


Suffonylurea receptor binrjing 




2 


17098 


0.00326252 


Calcium- and caimodulin-responsive adenyiate cyciase activity 


5 


8294 


0.00384257 


Positive reguiation of histone modification 




41 


31068 


0.00648655 


Negative reguiation of smali GlT*ase mediated 


signal transduction 


30 


51068 


0.00609217 


Mitochondrial respiratory chan 




5 


6746 


0.00636338 


Pyrimidine ribonucieoside binding 




4 


32661 


0.00722603 


Aging 




76 


7668 


0.00741073 


Carnitine transport 




4 


15879 


0.00883626 


Patiiway enriciiment 


Patliway name 


# of probesets 


Patiiway ID 


/^Value (attribute) 


Mucin lype 0-Giycan biosynthesis 


26 


Kegg pattiv^ay 


48 


0.00864224 


Ribosome biogenesis in eukaryotes 


75 


Kegg patiiway 1 00 


0.0276843 


Protein export 


19 


Kegg patiiway 22 


0.0642122 


Reguiation of autophagy 


28 


Kegg patiiway 9 


0.0966468 



Table 1 : Gene ontology and pathway analysis of Glut4 neurons. 

Top categories in each group of gene ontology analyses performed using RNA profiling data 

from flow-sorted Glut4 neurons. 



a less impressive increase of dopamine D2 receptor (Figure 1G). The 
expression of ion channels showed no particular pattern of enrichment 
in Glut4 neurons (Figure 1 J). We analyzed the expression of major sub- 
classes of ion channels (Figure 1 K— N), but they revealed no significant 
change. 

3.2. Characterization of InsR-deficient Glut4 neurons 
Our hypothesis was that Glut4 neurons are important to mediate CMS 
insulin signaling. To test this possibility, we interrogated the gene 
expression profiles of InsR-deficient Glut4 neurons using the same 
method described above (exon microarray hybridization of flow-sorted 
hypothalamic neurons), and performed electrophysiological studies. 
We used Tomato fluorescence to identify Glut4 neurons in WT and 
GIRKO mice bearing either two/one {i.e. Insi'°'^* or lnsr^'~) or no intact 
/nsfl allele (Insr^^^ or lnsr^'~). We collected Tomato-positive (i.e. Glut4 
neurons) and Tomato-negative (i.e. non-Glut4 cells) fractions from both 
WT and GIRKO mice (Figure 2A). Using DNA extracted post-sort from 
two independent sample sets (Figure 2B), we confirmed that Glut4 
neurons from GIRKO mice had minimum intact InsR allele. We 
compared gene expression profiles from sorted neurons (Figure 2C). 
Gratifyingly, gene set enrichment analysis of the array results showed 
that InsR was one of the most decreased mRNA (Supplemental 
Table 1), and that impaired expression of genes involved in the 
cellular responses to insulin topped the list of altered pathways, 
demonstrating the robustness of the method (Table 2). 
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Figure 1 : Glut4 neuron RNA profiling (A, J) Heatmap displaying the level of gene expression in differenf cafegories of neurotransmifter recepfors (A) and ion ctianneis (J) with control samples normalized to 1 , The iighter color indicates higher gene expression, the darker 
coior lower gene expression (n = 3), Individuai gene expression for each subcategory of neurotransmitter receptors (B-l) and ion channels (K-N) was expressed as bar graphs (n = 3) 
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Does InsR ablation in Glut4 neurons affect ttieir neurotransmitter 
profile? To answer ttiis question, we queried the dataset for major 
categories of neurotransmitter receptors. Interestingly, GIRKO Glut4 
neurons showed decreases in nearly all categories, as visualized by 
heatmaps (Figure 3A). We further analyzed the data by grouping major 
sub-classes of neurotransmitter receptors (Figure 3B-I). Even though 



most changes did not reach significance (Figure 3B-E, G-l), 
cholinergic receptors M4 and IV15 were significantly reduced in InsR- 
deficient Glut4 neurons (Figure 3F). Of note, the predominant 
GABAergic receptor tended to be increased (Figure 3H). Expression of 
ion channels was reduced to different extents in InsR-deficient Glut4 
neurons (Figure 3J). While most calcium, chloride, and sodium 
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Gene set description 


# of marlters 


p-Vaiue 


Description 


Cellular response to peptide hormone stimulus 


91 


0 


KO down vs WT 


Response to Insulin stimulus 


97 


0 


KO down ys WT 


Anion homeostasis 


36 


0.0185 


KO up vs WT 


Rho guanyi-nucieotide exchange factor activity 


52 


0.020 


KO down ys WT 


Regulation ot transcription Initiation 


15 


0.020 


KO down ys WT 


Epldemiai growtii factor receptor binding 


20 


0.0208 


KO down ys WT 


Celiular protein complex assembly 


152 


0.0213 


KO down ys WT 


Regulation ot cAMP metabolic process 


95 


0.0227 


KO down ys WT 



Table 2: Gene ontology analysis of RNA profiling data from WT and KO Glut4 neurons. 
Gene enrichment analysis showing the most significantly altered cellular pathways using RNA 
profiling data from flow-sorted WT and InsR-KO Glut4 neurons. 



channels showed comparable or slightly reduced expression In InsR- 
deficient Glut4 neurons (Figure 3K-IVI), K-channel KcnjIO was 
significantly reduced in InsR-deficient Glut4 neurons (Figure 3N). 
Other top pathways with decreased expression In Glut4 neurons of 
GIRKO mice included guanyl-nucleotlde exchange factor (GEF) activity, 
regulation of transcription initiation, cellular protein complex assembly, 
and regulation of cAWIP metabolic processes. Conversely, anion ho- 
meostasis was the most significantly up-regulated gene set In GIRKO 
Glut4 neurons (Table 2). 

3.3. Altered firing properties of InsR-deficient Glut4 neurons 
The changes In anion homeostasis and neurotransmltter/ion channel 
expression would be expected to result in altered firing properties of Glut4 
neurons. To test this hypothesis, we performed electrophysiological ex- 
periments (Figure 4A-D). We found that patch-clamped InsR-deflclent 
Glut4 neurons had significantly reduced firing frequency (Figure 4C and 
D), despite similar resting membrane potential, and action potential 
amplitude compared with WT neurons (Figure 4A and B). 

4. DISCUSSION 

The l<ey findings of our work are that hypothalamic Glut4 neurons are 
characterized by their enrichment in genes required for olfacto-sensory 
and metabolic processes, as well as the muscarinic cholinergic re- 
ceptor M2. Moreover, ablation of InsR results in loss of insulin signal 
transduction properties, and decreased KAjp-channel activity, as 
revealed also in patch-clamp experiments. The data are consistent 
with a primary metabolic coupling function of these neurons and their 
insulin responsiveness. 

Although we found general enrichment of neurotransmitter receptors 
and ion channels in Glut4 neurons, we did not observe enrichment of 
particular sub-classes. This indicates that hypothalamic Glut4 neurons 
are a heterogeneous group of neurons and lil<ely have multiple 
metabolic functions. Indeed, the distribution of Glut4 neurons to 
multiple neuroanatomical sites in the brain, including cortex, hippo- 
campus, and hypothalamus [7] suggests that they have diverse bio- 
logical functions. 

The findings are largely consistent with the hypothesis that the pres- 
ence of Glut4 in these neurons is a hallmarl< of metabolic functions. We 
should emphasize that in our studies we didn't probe the function of 
Glut4 itself in this process, let alone whether it responds to insulin in 
ways similar to muscle and fat Glut4. The present findings begin to 
explain why ablation of InsR in Glut4-expressing tissues has much 
more profound consequences on glucose homeostasis than ablation of 
InsR limited to "canonical" insulin-responsive tissues (i.e. muscle and 
fat) [7,11]. 

The Identification of olfacto-sensory pathways in Glut4 neurons sug- 
gests that one potential function of this neuronal population is to 



connect sensorial and metabolic aspects of energy balance. These 
findings are reminiscent of findings using retrograde mapping of leptin 
receptors in the hypothalamus, demonstrating that they project to 
neurons in the premammilary tract that receive olfactory inputs [12]. 
IVIoreover, intranasal delivery of insulin is l<nown to affect olfactory 
behavior via Kv1.3 channels [13]. In this regard, it's worth noting that 
six of the top ten genes downregulated in InsR-deficient Glut4 neurons 
encode olfactory receptors (Supplemental Table 1). The Quebec Family 
Study has found a genetic association between a cluster of olfactory 
receptors on human chromosome 19 and eating behavior [14]. Our 
data provide a potential lead on how olfacto-sensory cues might 
become integrated in energy balance and peripheral metabolism. 
While Glut4 neurons don't appear to be enriched in a specific subset of 
neurotransmitters, they show significantly increased levels of the M2 
muscarinic acetylcholine receptor. This receptor subtype in the hip- 
pocampus has been Implicated in short- and long-term synaptic 
potentiation as well as neuronal plasticity that affect behavior [15], but 
considerably less is l<nown about Its role In the hypothalamus [1 6]. Our 
data provide a rationale to interrogate its metabolic contributions more 
specifically in this anatomical location. 

Our gene enrichment analysis indicates the presence in Glut4 neurons 
of a gene cluster commonly associated with adipocyte differentiation 
that includes Ppar-y, Tcf7l2, Klf4, C/ebp-a and -fi, LIpin, and several 
others. Interestingly, Ppar-y activation in the CNS has been linl<ed to 
decreased insulin sensitivity [17]. Similarly Tcf7l2, a diabetes sus- 
ceptibility gene [1 8] that had been originally linl<ed to fi-cell function, is 
expressed at high levels In the midbrain [19]- consistent with the 
present data - and has been found to regulate glucagon gene 
expression [20]. 

The transcriptome analysis of InsR-deficient hypothalamic Glut4 neu- 
rons supports key conclusions of our study and is consistent with an 
important role of these neurons in CNS Insulin action. We have listed 
the top 1 00 most down- or up-regulated genes In InsR-deficient hy- 
pothalamic Glut4 neurons here (Supplemental Tables 1 and 2). First, 
gene ontology analyses uniformly Indicate a loss of insulin and peptide 
hormone signaling pathways, as well as cAMP and transcriptional 
responses. The only class of genes whose expression increases 
following InsR ablation is involved in anion homeostasis. Interestingly, 
increased plasma anion concentrations are commonly seen In insulin- 
deficient states (e.g., diabetic ketoacidosis), and are thought to reflect 
increased generation of ketone bodies as an alternative energy source 
for neurons, possibly via glial cells — although this hypothesis remains 
controversial [21]. The present findings are consistent with the idea 
that, In conditions of impaired insulin signaling, there is an Increase in 
the neuron's ability to utilize anions, possibly as a survival mechanism 
when glucose utilization is impaired. 

As mentioned above, several olfactory receptors are sensitive to InsR 
ablation. Interestingly, one the top downregulated genes in InsR- 
deficient Glut4 neurons encodes a ribosomal protein (rsg24), whose 
decreased expression In the CNS has been linked to type 2 diabetes 
[22]. The Crh receptor also made the list of top downregulated genes, 
suggesting altered CNS stress response with Insulin resistance. 

5. CONCLUSION 

In sum, our data define hypothalamic Glut4 neurons as a cell popu- 
lation with characteristics of hormone responsiveness, and with the 
potential to link Insulin signaling to olfacto-sensory signaling. By 
systematic analyses of Glut4 neurons in additional CNS areas, it should 
be possible to gradually piece together an Integrate anatomic and 
functional map of insulin action in the CNS. 
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Figure 3: Glut4 neuron RNA profiling from WT and KO Giuf4 neurons (A, J) Heatmap dispiaying fhe ievei of gene expression in different cafegories of neurotransmitfer receptors (A) and ion channels (J) with control samples normalized to 1 . The lighter color indicates 
higher gene expression, fhe darker color lower gene expression (n = 3). Individual gene expression for each subcategory of neurotransmitter receptors (B— I) and ion channels (K-N) was expressed as bar graphs [n = 3). 
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Figure 4: Eiectrophysiological characterization ot Glut4 neurons from WT and KO mice (A— D) Electrophysioiogicai properties of VIvlH Giut4 neurons. Patch-ciamp recordings of the spontaneous firings 
of VMH Glut4 neurons in WT and GIRKO mice, (A) Resting membrane potentiai (RMP). (B) Ampiitude of action potential. (C) Firing frequency. (D) Representative traces. Data show means ± SEM 
(n= 19-24) r = p< 0.01), 
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